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Abstract 

Thermal contact conductance test data at high 
vacuum were obtained from two Armco iron specimens 
having smooth, large, radii of curvature, convex, 
one -half wave length surfaces. The data are com- 
pared with calculations based on two macroscopic 
elastic deformation theories and an empirical ex- 
pression. Major disagreement with the theories and 
fair agreement with the empirical expression result- 
ed. Elastic deformation of all the contacting sur- 
faces was verified from surface analyzer statistics. 
These results indicate that the theoretical assump-, 
tion of macroscopic elastic deformation is inade- 
quate for accurate prediction of heat transfer with 
light loads for Annco iron specimens similar to 
those used in this investigation. 

Introduction 


Smooth metal surfaces having large radii of 
curvature joined together under light load by riv- 
ets, bolts, etc. are commonly found in spacecraft. 
Such joints act as paths through which heat is 
transferred. The theories of A. M. Clausing and 
B. T. Chaol 1 ) and B. B. Mikic and W. M. Rohsenow^ 2 ) 
are frequently used to predict the coefficient of 
thermal contact conductance for contacting surfaces 
under high vacuum conditions at light loads. Two 
of the assumptions made in these theories may not 
be strictly applicable to specimens having smooth 
contacting surfaces with large radii of curvature. 
The first assumption specifies that macroscopic 
elastic deformation occurs during the formation of 
such joints with an associated effect on the heat 
transfer rate. The second specifies that the mi- 
nute contact spots which form during the formation 
of a joint are uniformly distributed over the con- 
tacting area. Discussions in Refs. 1 and 2 con- 
cerning test results indicate possible invalidity 
of these assumptions. Clausing and Chaov^-/ indi- 
cate consistently poor theoretical agreement with 
data from several specimens having smooth large 
radii of curvature. Also, they note that plastic 
deformation occurred for three magnesium specimens 
which had smooth surfaces with large radii of cur- 
vature. In Ref. 2 Mikic and Rohsenow indicate that 
their theoretical expression might have serious 
limitations because it is dependent on the uniform 
distribution of the contact spots over the contact- 
ing, area. Mikic and Rohsenow conclude that if the 
distribution of the contact spots were not uniform 
it would have the same effect on the conductance as 
an equivalent type of waviness. They present ex- 

f arimental data from several rough nominally flat 
i.e., no waviness) specimens to substantiate the 
effect. Because of the lack of agreement between 
experimental data and the theories of Ref. 1 and 
other similar theories, L. S. Fletcher and D. A. 
Gyorog(3) developed an empirical expression for the 
coefficient of thermal contact conductance. They 
show correlation between their expression and ex- 


perimental data and the data of seven other inves- 
tigators including Clausing and ChaoC 1 ) to within a 
mean deviation of 24 percent. Because of this cor- 
relation with an immense amount of data, their work 
is of particular interest here. 

This paper presents experimental test data ob- 
tained from two Armco iron cylindrical specimens 
having smooth, convex, one -half wave length contact- 
ing surfaces with large radii of curvature which 
show the variation of the coefficient of thermal 
contact conductance with apparent contact pressure 
in the light to moderate range (0.45x10® to 3.44>d_0® 
N/m 2 ) . The data were gathered during two loading 
and unloading sequences of each specimen at a con- 
stant interface temperature (369 K) under high vac- 
uum conditions (l.2xl0 - ? torr). Results are com- 
pared with the theories of A. M. Clausing and B. T. 
Chao(l) and B. B. Mikic and W. M. Rohsenow. C 2 ) 

Also, a. comparison is made with the empirical solu- 
tion of L. S. Fletcher and D. A. Gyorog.(3) 

Symbols 

radius of macroscopic circular contact 
area 

apparent contact area 

radius of cylindrical specimen 

half-wave height of specimen' s spheri- 
cal cap 

surface parameter 

modulus of elasticity 

flatness deviation 

function representing the effect of 
neighboring contact points on the 
flow of heat through a point con- 
tact 

H micro-hardness of the softer of the two 

contacting surfaces 

h coefficient of thermal contact conduc- 

tance 

k harmonic mean of the coefficient of 

thermal conductivity 

1 length of upper or lower portion of a 

specimen 

L s one-half the spherical surface wave 

length 

n number of contacts spots per unit area 


a 

A 

b 

d 

d' 

E 

FD 

f(x L ) 
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a 

P 
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apparent contact pressure 

heat transfer rate 

radius of curvature of the contact 
surface 

rms roughness 

temperature at interface of contact 

applied load 

the constriction ratio 

contact area ratio 

coefficient of thermal expansion 

effective surface parameter represent- 
ing the gap between the contacting 
sur faces 


(g^ax) maximum average roughness height ob- 

tained from total profile traces 

2 

€ ratio of real to apparent contact area 


*eff 


effective diameter of the contour area 


v 


Poisson' s ratio 


*(*eff) 


functions which represent approximate 
finite series solutions of the tem- 
perature field at the contact plane 


Subscripts: 


resistance is defined as the reciprical of the prod- 
uct of the apparent contact area. A, end the coeffi- 
cient of thermal contact conductance, h. The appar- 
ent contact area is equal to the cross-sectional 
area of their cylindrical model. The theory applies 
to contacting surfaces loaded in the light to mod- 
erate range at high vacuum conditions, where it is 
assumed that elastic deformation of the macroscopic 
area of contact of each member occurs. Their model 
is shown in Fig. 1 which is similar to one appear- 
ing in Ref. 1. The theory considers the contact 
region between two cylinders of length L and 
identical radius b placed end to end. The ends 
of the cylinders in the contact plane have radii of 
curvature r^_ and r g (top and bottom, respect- 
ively) as shown in Fig. 1. A single large, ’’macro- 
scopic, circular contact area of radius a is 
assumed to be concentric with that of the apparent 
contact area. Thus the noncontact region is an 
annulus whose inside and outside radii are a and 
b, respectively. Inside the macroscopic contact a 
high density of microcontact areas is assumed. The 
microcontact areas are assumed circular, of the 
same radius, and uniformly distributed over the 
single macroscopic spot. To determine the change 
of macroscopic contact area with load, the surfaces 
are assumed to approximate spherical caps; thus the 
Hertzean solution! 4 ) can be applied. Therefore, 
the Clausing and Chao macroscopic theory predicts 
the coefficient of thermal contact conductance be- 
tween surfaces having a spherical contour of one- 
half wave length. The form of the solution taken 
from Ref. 1 and used in this paper is 


hb = 

k «f(x L ) 


(1) 


1 top portion of specimen 


where the limits of its application are 


2 bottom portion of specimen 

i interface 



no load 


Theoretical and Empirical Analyses 

The theories of Clausing and Chao(^) and Mikic 
and Rohsenow( 2 ) both consider distinct contribu- 
tions to the thermal conductance resulting from the 
plastic and elastic deformation of the contacting 
surfaces. Specifically, the roughness of a sur- 
face is assumed to deform plastically while ■ the 
waviness is assumed to deform elastically. Asper- 
ities produced by machine tools in the forming of 
surfaces are called roughness (microscopic irregu- 
larities); and surface irregularities having wave 
lengths generally greater than 0.080 centimeters 
axe called waviness (macroscopic irregularities). 
Though both theories consider the effects of these 
deformations the final form of the Clausing and 
Chao solution expresses the roughness contribution 
only approximately, therefore, only the waviness 
contribution of their theory (macroscopic theory) 
will be described here. 


Clausing and Chao's Macroscopic Theory 

In 1963 A. M. Clausing and B. T. Chao(^) 
presented their macroscopic theory of 
constrictive resistance. The constrictive 


k is the harmonic mean of the coefficient of 
thermal conductivity of the two contacting materials 
evaluated at their contact interface, and f(x,) 
represents the effect of neighboring contact points 
on the flow of heat through a point contact. f(x L ) 
is evaluated from 

f(^) = 1 - 1.40925 ^ + 0.29591 x^ + 0.05254 

+ 0.02105 x^ + 0.01107 xj? 

where x^ is the constriction ratio or equivalently 
the square root of the ratio of the contact area to 
the apparent contact area. For two spherical sur- 
faces in contact x^ is. evaluated from 



where v-|_ and v g axe Poisson' s ratio which for 
most metals one may assume vf = v| = 0.1, E-j_ 
and E g are the moduli of elasticity, W is the 
applied load, and r^ and r g are the , radii of 
curvature of the spherical surfaces. In the deri- 
vation of this equation Clausing and Chao indicate 
that the dimensions of the bodies in contact were 
assumed large in comparison with the radius of the 
macroscopic contact area. For this reason, Eq. (2) 
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can be applied to the macroscopic contact model 
only for small values of x^. On the basis of 
several tests reported by them, Eg,. (2) was found 
applicable for xr less than approximately 0.65. 

The result of a study performed by them to deter- 
mine the influence of the length of each portion of 
the cylindrical specimen on the macroscopic conduc- 
tance indicated that for values of L/b greater 
than approximately 0.6 the conductance was unaf- 
fected. 

Mikic and Roshsenow' s Theory 

In 1966 B. B. Mikic and W. M. Rohsenow^ 2 ) pre- 
sented methods for predicting the coefficient of 
thermal contact resistance for four different types 
of surfaces, however, only the rough and spheri- 
cally wavy type will be considered here. The ther- 
mal contact resistance is presented in Ref. 2 as 
the reciprical value of the thermal contact conduc- 
tance. Mikic and Rohsenow consider the surfaces to 
be pressed together with intimate contact occurring 
only at a discrete number of locations called con- 
tact spots which are uniformly distributed over the 
contact surface. The spots are all assumed to be 
circular, flat, and have the same area. These spots 
introduce the concept of surface asperity roughness 
into their model. The contact spots are assumed 
confined inside contour areas which in turn have 
their own distribution. The contour areas intro- 
duce the concept of waviness into the model. For 
the specific case considered here, the contour 
areas are assumed to approximate spherical segments. 
The surface analysis assumes the existence of an 
ensemble of the surface roughness -only profiles, 
all taken from one surface, from which one can de- 
duce statistical properties of the surface and, 
therefore, a random process which possesses a prob- 
ability density function. The random process is 
assumed stationary, i.e., the statistical properties 
of the ensemble of surface profiles are invariant 
under arbitrary displacement over the surface. In 
addition, the probability density of the height and 
slope of the profiles are independent and the sur- 
face height is assumed to have a standard normal 
(Gaussia n) distribution. From this analysis the 
number of contacts per unit area and hence the real 
area of contact and the effective diameter of the 
contour areas are determined. In the thermal anal- 
ysis the temperature distribution and implicitly 
the thermal contact resistance is specified by the 
Laplace differential equation (for steady state 
conditions and thermal conductivity independent of 
temperature). Because heat is assumed to flow only 
through flat contact spots, mixed boundary condi- 
tions exist at the surface and the temperature field 
at the contact plane is obtained by the method of 
superposition of an infinite number of heat sources 
equally spaced on the contact surface which extends 
to infinity. 

The form of the prediction of the conductance 
taken from Ref. 2 for two rough and spherically 
wavy surfaces exposed to a vacuum environment is 

jfej ~ ^(^eff^s 

ke-v/’W kX e ff 


tively. In Eq. (3) h is the coefficient of ther- 
mal contact conductance, € 2 is the ratio of the 
real to apparent contact area, n is the number of 
contact spots per unit area, L s is equal to one- 
half the spherical surface wave length, is the 

effective diameter of the contour area, and k is 
the harmonic mean of the coefficient of thermal 
conductivity of the two contacting materials evalu- 
ated at their interface of contact. 4>( e ) and 
(j>(k e ff) are factors which represent the approxi- 
mate finite series solution of the temperature field 
at the contact plane. For a more detailed descrip- 
tion of the parameter which involves statis- 

tical quantities, refer to Ref. 2. The effective 
to apparent area ratio of a contour area, c 2 , is re- 
lated to the apparent contact pressure by assuming 
that the surface asperities deform plastically. 
Therefore, 

H(Aeff) (4) 

where P is the apparent contact pressure and H 
is the micro hardness of the softer of the two con- 
tacting surfaces. 

In a discussion of the limitations of the pre- 
diction of the conductance, the authors point out 
that the assumption of the same size and uniformity 
of distribution of the contact spots over the appar- 
ent contact area might have serious limitations on 
the accurate prediction of the coefficient of ther- 
mal contact conductance especially at light pres- 
sures. The degree of nonuniformity has the same 
effect on the conductance as some equivalent type 
of waviness. Experimental results from several 
nominally flat specimens were presented to demon- 
strate this effect. 

Fletcher and Gyorog's Empirical Solution 

In 1970 L. S. Fletcher and D. A. Gyorog( 3 ) pre- 
sented an expression for the prediction of the ther- 
mal conductance of similar contacting metal surfaces 
in a vacuum environment. Fletcher and Gyorog point 
out that although theoretical and experimental 
studies concerning heat transfer phenomenon in solid 
materials have resulted in many useful techniques 
for analysis, few attempts have been directed at 
the correlation of existing experimental data and 
thus the prediction of contact conductance. In 
their study the experimental results for smooth, 
medium, rougu, and smooth to rough contacts of 
aluminum, stainless steel, brass, and magnesium 
were used to provide information for the develop- 
ment of an empirical solution. 


To classify the contact surfaces used in their 
study, a surface parameter, d' , was defined in terms 
of its rms roughness and flatness deviation. The 
contact between surfaces was assumed to occur at 
the mean surface height of the smoother surface. 

The contact surface parameter is defined as 


d- H (FD + 2RD) rough - |(FD + 2RD) smooth 
surface surface 


( 5 ) 


where FD is the flatness deviation, and RD is 
the rms roughness. 


The left and right hand terms in the denominator of 
Eq, (3) represent the effect due to the plastic 
deformation of the surface roughness and the elas- 
tic deformation of the surface waviness, respec- 


From previous work by Fletcher a theoretical 
expression for the conductance, h, was developed 
for a single contact and is of the form 


3 



2L_ (6) 

. .Urn 1 - « 

where a is the contact area ratio, 5 is an effec- 
tive surface parameter representing the gap between 
the contacting surfaces, and is the mean value 

of the thermal conductivity. ' From- Eq. (6) and 'an 
analysis of their experimental" data in dimension- 
less form, the value of a no-load surface parameter 
6 0 for each of the surface configurations was de- 
termined. The no-load surface parameter & 0 re- 
sulted by assuming that the differences in the test 
data were due only to differing no-load surface con- 
ditions. The values of 5 0 for all specimens were 
plotted as a function of the contact surface param- 
eter, d' . From this it was evident that the no- 
load surface parameter. So, was related to the con- 
tact surface parameter, d' , and thus to the meas- 
ured quantities of flatness deviation and roughness. 


This work along with some further analysis of 
the no-load surface parameter, 5 0 , and the effec- 
tive surface parameter, 6, led to the final form of 
the conductance expression 


^m 


177(PpT m b)/(E&, 


o) 5. 


22x10 ' 6 

b 


+ 0.036 


PpT, 


10.56 ' 

T 


(?) 


where P is the apparent contact pressure, p is 
the coefficient of thermal expansion, T m is the 
mean temperature at the interface of contact, b is 
the radius of the specimen, and E is the modulus 
of elasticity. 


The suggested form of 8 0 is given by 


6 = 20,45 + 8.06Xl0~ 2 d' - 1.58x10 _5 (d' ) 2 
o 

+ 1.36xl0~ 9 (d’) 3 


where 6 0 and d.' , the contact surface parameter, 
are in microinches. 

Fletcher and Gyorog note that the expression 
for the conductance, i.e,, Eq. (7), includes pos- 
sible radiation effects as the apparent contact 
pressure is reduced to zero. Also they point out 
that the expression correlates the experimental 
data of their investigation and those of seven other 
investigators whose data comprised approximately 
400 points. These data include a range of mean 
interface temperatures varying from 117 K to 532 K, 
apparent contact pressures of 6.89X10 4 to 4. 82x10 7 
N/m^, surface flatness deviations of 38.1 to 11420 
microcentimeters, and surface roughness of 7.6 to 
305 microcentimeters. Data for aluminum, stainless 
steel, brass, magnesium are included. The pre- 
diction approximates the data within a mean devia- 
tion of 24 percent or less. 

The following is a brief discussion of the 
applicability of the empirical expression developed 
by Fletcher and Gyorog to the Armco iron test re- 
sults discussed herein. The test specimens' sur- 
face characteristics obtained by Fletcher and 
Gyorog are listed in Ref. 3, table I. These char- 
acteristics are indicated as the average flatness 
deviation and the average roughness. In all cases 


but one, the average flatness deviation is 5 times 
larger than the average roughness of each surface. 
In.the single exception the average flatness devia- 
tion is greater than twice the average roughness. 
This result indicates that none of the surfaces were 
nominally flat (i.e., not wavy) and at least two 
thirds of them can be considered smooth and wavy. 
Therefore, the expression derived by Fletcher and 
Gyorog would be expected, to apply to surfaces which 
are characterized as wavy with some degree of rough - 
"ness, but would not apply to surfaces which are 
characterized as being nominally flat. Their ex- 
pression, therefpre, is particularly applicable to 
the Armco iron test data presented herein. 

Experimental Apparatus and Procedure 

Facility . - The two Armco iron specimens were 
tested in the facility shown in Fig. 2. The tests 
were performed in a bell jar containing a column 
composed of a tungsten resistance heated cylindrical 
heat source and its supporting structure, the test 
specimen, and a water cooled heat sink. The lower 
portion of each specimen was especially designed, 
instrumented, and calibrated against a 99.993 per- 
cent pure aluminum (by assay) bar to serve as a 
heat meter. A pheumatic loading cylinder was locat- 
ed on top of the bell jar and was used to vary the 
pressure at the contact interface of the test speci- 
men. The heat generated by the heater was trans- 
ferred across the test interface, through the heat 
meter, and finally to the water cooled heat sink. 

The bell jar operates in the 10 - 6 torr pressure 
range. 

Specimens . - The Armco iron cylindrical speci- 
mens were 2.54 centimeters in diameter. Their upper 
portions were 6.97 centimeters long and their lower 
portions were 12.7 centimeters long. 

Surface finish . - The specimens' contacting 
interface surfaces were produced by lapping and had 
no lay. Table I presents the average values of the 
specimens maximum roughness profile heights, 

( &max) ave (defined in ref. 5) and their arithmetic 
average heights A. A. (defined in American Standard 
ASA B 46.1 - 1962 ASME). The contacting surfaces 
of the specimens were found to approximate sections 
of spherical caps having one -half wave length, 
therefore, the average crest to trough wave height, 
d, (based on radius of curvature) for the surfaces 
is also shown. 

Thermometry . - The thermocouples were located 
in the specimens as shown in Fig. 2 and read out on 
a potentiometer. Type 30 gage, special error lim- 
its, chromel-alumel thermocouple wire, butt welded 
at its ends, was used. These thermocouples have an 
accuracy of ±1.1 K as guaranteed by the manufacture 
in the temperature range of their use here. After 
the tests were completed these thermocouples were 
calibrated in place on the specimens and were found 
to be within the manufacture's tolerance for the 
range of temperatures -occurring during the tests. 

The couples were peened into slots milled into the 
circumference of each specimen. The slots were 
approximately 0.076 centimeter long, 0.033 centi- 
meter wide, and at their deepest point were approxi- 
mately 0.051 centimeter deep. The thermocouple's 
insulated lead wires were wrapped around the circum- 
ference of the specimens to minimize conduction 
errors in the thermocouple wire. Four couples were 
located in the upper portions of the specimens at 
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1.27, 2.54, 3.81, and 5.08 centimeters from the 
interface of contact. Six couples were .located in 
the specimens' lower portions at 1.27, 2.54, 3.81, 
6.35, 8.89, and 11.5 centimeters from the interface 
of contact. 

Test Procedure 

The contact conductance tests were performed 
by varying the interface contact pressures of the 
spec ime ns in a monotonic increasing then decreasing 
manner through two complete cycles from approxi- j 
mately 4.5x10° to 3.5x10° N/m 2 (65 to 500 psi) 
while maintaining their mean interface temperatures' 
nearly constant at approximately 367 K. 

To gain an idea of the time required to estab- 
lish equilibrium between changes in pressure and/or 
temperature the test specimens were exposed to 
steady-state test conditions for periods ranging 
from 50 to 180 hours between each point while moni- 
toring the variation of the conductance with time. 
Results from this procedure indicated a minimum of 
50 hours was required to reach a stable condition 
between pressure changes. 

Before placing the test specimens in the test 
column of the bell jar they were cleaned with ace- 
tone and baked out in a separate vacuum facility at 
10"® torr, through eight 30 minute heating and 
cooling cycles. Each cycle varied, from room tem- 
perature to approximately 450 K, well below the re- 
crystallization temperature of Armco iron. 

Data Reduction . - The coefficient of thermal 
contact conductance of each specimen, h, was experi- 
mentally determined from the • relation 

i=-h(t ij 2 -t i5l) (e) 

where Q/A represents the heat transfer rate per 
unit area across the interface of contact and was 
determined by a heat balance performed across each 
specimen' s heat meter, discussed in the test appa- 
ratus section. The interface temperature differ- 
ence (ti 2 - ti i) was determined by the linear 
extrapolation of uhe thermocouple data in the upper 
and lower portions of each specimen to the inter- 
face of contact. 

A simplification in this procedure was made by 
determining the thermal gradient in each portion of 
a specimen from two thermocouple readings in the 
same manner as that used by Fried. ( 6 ) The gradi- 
ents determined by this technique were continually 
checked by plots including all of the thermocouple 
data from each specimen. The error in temperature 
difference across the interface introduced by this 
method was less than 1.1 K from a best line (drawn) 
through all the thermocouple data from each speci- 
men. 

Theoretical Determination of h 

Before the theoretical predictions of Refs. 1 
and 2 and the empirical expression of Ref. 3 could 
be applied to determine the coefficient of thermal 
contact conductance for the test specimens consid- 
ered here, it is necessary to obtain the roughness 
profile, total profile, and waviness of each of the 
contacting surfaces. These surface texture quanti- 
ties were determined by using a Brush surfanalyzer 


1200 system. A description of this system, its 
operation, and use is given in Ref. 7. 

The surface texture quantities, obtained di- 
rectly from the surfanalyzer traces, were used in 
the theoretical prediction of Clausing and ChaoW 
and also the empirical solution of Fletcher and 
Gyorog. ( 3 ) However, because of the nature of the 
theoretical prediction of Mikic and Rohsenow' 2 ' it 
was necessary to use a digital computer program for 
the case of spherically wavy, rough surfaces in 
contact. A description of the computer program and 
its use is given in Ref. 7. 

The material properties appearing in the theo- 
retical and empirical expressions used in the cal- 
culations were obtained from values in the litera- 
ture with the exception of the coefficient of ther- 
mal conductivity, k, and the microhardness, H. The 
variation of k with temperature was determined 
from calibration data taken in the test facility 
under high vacuum conditions for both Armco iron 
specimens using a 99.993 percent (by assay) pure 
aluminum heat meter. The results of this calibra- 
tion indicate agreement with data from the litera- 
ture^) to within 1.5 percent. The contact sur- 
faces' microhardness, H, was measured using a 
Reichert microhardness tester. These values were 
found to be larger by approximately 1.7 times the 
Meyers hardness values appearing in the literature. 

Results and Discussion 

Experimental cyclic test data are presented 
for two Armco iron cylindrical specimens having 
smooth, convex, spherical, one-half wave length 
contacting surfaces which show the variation of the 
coefficient of thermal contact conductance with 
apparent joint contact pressure in the light to 
moderate range (0.45x10° to 3.44XL0 6 N/m 2 ). These 
data are compared with the theoretical predictions 
of Clausing and 0hao(-'-) and Mikic and Rohsencw. ' 2 ' 

In addition, the data are compared with the empiri- 
cal solution of Fletcher and Gyorog.'^' 

Cyclic Test Data 

Figure 3 presents the experimental test data 
obtained from specimen 1. This specimen was loaded 
through two cycles. Its mean interface temperature 
was kept constant at 369 K with one exception noted 
at the end of the first loading where it was raised 
at constant pressure (350 N/m 2 ) to 477 K. The 
interface remained at this temperature for 50 hours. 
Then while maintaining constant pressure it was 
lowered to 369 K and the pressure cycling was con- 
tinued. The data obtained for the first loading 
followed the generally recognized trend exhibited 
when macroscopic elastic deformation of contacting 
surfaces takes place. All data following the first 
loading fell reproduciably on a curve having a 
slope slightly greater than one which is character- 
istic after moderate loading pressures have been 
applied. For the smallest apparent contact pres- 
sures, the measured thermal conductances were con- 
sistently lower than those obtained initially. 
Therefore, it may be surmised that the specimen 
took a permanent set after the initial loading. 

The results from specimen 2 obtained under 
almost identical experimental conditions are shown 
in Fig. 4. These data appear in Ref. 5 in a 
slightly different form. Unlike specimen 1, how- 
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ever, specimen 2 was not exposed to a temperature 
excursion. The data of specimen 2 follow the ac- 
cepted trend for elastic deformation. However, 
closer examination at the lowest loading pressure 
shows a consistent and reproducible decrease in the 
conductance perhaps indicating an increasing per- 
manent set. 

The contacting surfaces were exam i ned before 
and after each test with a Brush surfanalyzer. 

Table I presents each specimen' s surface half-wave 
height, d, obtained before and after the tests. 

The half-wave height of the top portion of speci- 
men 1 increased bv 335x10"® centimeters and the 
bottom by 210>d.0"° centimeters. The top of speci- 
men 2 increased by 149x10 "® centimeters and the 
bottom by 122x10"° centimeters. Thus both speci- 
mens deformed plastically with the attendent 
changes in the coefficient of thermal contact con- 
ductance noted in the test data. The interface 
temperature excursion apparently caused specimen 1 
to deform more than specimen 2. 

In Ref. 1 Clausing and Chao obtained test re- 
sults from a total of 22 specimens fabricated from 
brass, magnesium, stainless steel, and aluminum. 

The data obtained from one of three magnesium spec- 
imens indicate that its total equivalent flatness 
deviation increased from before (89x10"® cm) to 
after the test (203.2x10"® cm). Where the total 
equivalent flatness deviation, referred to in 
Ref. 1, is equal to the sum of each portion of a 
specimen' s half-wave surface height, as used in this 
paper. With the exception of the three magnesium 
specimens, no permanent set was detected among 
these specimens. This was determined by analyzing 
each contact surface before and after a test using 
optical interference fringe techniques. 

Comparison of Data With Theory 

Figures 5 and 6 show the curves determined by 
applying the theories of Refs. 1 and 2 and the 
empirical expression of Ref. 3 to the test data 
obtained during the first loading of specimens 1 
and 2, respectively. The calculations were made 
using the surface statistics listed in table I 
obtained before the tests were conducted. The co- 
efficients determined from both theories are ap- 
proximately 3 and 1.5 factors larger than the test 
data obtained from specimens 1 and 2, respectively. 
The values determined from the empirical expres- 
sion are approximately 35 percent lower than the 
test data obtained from the specimens. The dis- 
agreement between specimen' s 1 and 2 test results 
and the theories of Refs. 1 and 2 are similar to 
several cases discussed in Ref. 1. Clausing and 
Chao noted better agreement occurred for their 
specimens having rough surfaces and a total equiva- 
lent flatness deviation which varied from 
304.8x10"® to 2413x10"® centimeters than for sur- 
faces having a total equivalent flatness deviation 
which varied from 38.1x10"® to 101. 6x10"® centi- 
meters (the range of the Armco iron data presented 
here). Disagreement with theory occurred for speci- 
mens made from magnesium, stainless steel, and 
aluminum, but not for the brass specimens which had 
the largest total equivalent flatness deviations. 

As an explanation for the disagreement. 

Clausing and Chao suggest that visible oxide films 
were inhibiting the. flow of heat across their 
specimens' contacting surfaces. Wo such oxide 


films were visible on the Armco iron specimens test- 
ed. here, however. A discussion by Mikic and 
Rohsenow in Ref. 2 may serve as a second explanation 
for the large difference between theory and experi- 
ment noted here and in Ref. 1. In the theoretical 
and empirical analyses section of this paper it is 
indicated that the theoretical expression of Mikic 
and Rohsenow, as mentioned by them in Ref. 2, might 
have serious li m itations on the accurate prediction 
of the conductance. This is dependent on the degree 
of nonuniformity of the distribution of the contact 
spots over the apparent contact area. They conclude 
that if the distribution of the contact spots were 
nonuniform it would have the same effect on the 
conductance as an equivalent type of waviness. 

Though their discussion was in regard to rough nomi- 
nally flat surfaces in contact, the remarks could 
also apply to the smooth, wavy case where the dis- 
tribution of contact spots over the contour area is 
nonuniform. Like Ref. 2 it is assumed in Ref. 1 
that microscopic contacts are distributed uniformly 
over the contacting surfaces. Referring to the 
surface statistics obtained from specimens 1 and 2, 
table I shows that each contact surface had a sma.1 1 
number of very tall asperities (Bmax)ave (i.e., 
when compared with each surface' s arithmetic aver- 
age asperity height). As a result, it is probable 
that their distribution was not uniform over the 
contact surface, particularly at light loads. Thus 
as discussed by Mikic and Rohsenow and demonstrated 
by the test data of Clausing and Chao and the tests 
conducted here, the assumption that the contact 
spots are uniformly distributed over the contacting 
surfaces, made in both Refs. 1 and 2, may be too 
restrictive for smooth, convex, large radii of 
curvature surfaces i.e., surfaces having small half- 
wave heights. 

It should be noted here that the test results 
of specimen 2 were compared in Ref, 5 with the 
macroscopic theory of Clausing and Chao using the 
surfanalyzer statistics obtained before the test 
was performed. However, a mathematical error was 
made in the application of the theory, therefore, 
the theoretical curve presented in Ref. 5 is in- 
correct. The curve presented in Fig. 6 herein and 
labeled Clausing and Chao, macroscopic theory(-*-) is 
the correct curve. 

Figures 7 and 8 show the results of applying 
the theories and the empirical expression based on 
the surface statistics obtained from the specimens 
after the tests were made. Though the curves de- 
termined from the theories are in good agreement 
with the data for both specimens, the trend of the 
curve, in the case of specimen 1, is in poor agree- 
ment with the trend of the data. Also the use of 
the surface statistics obtained after the tests is 
not permissible since plastic set has been verified 
by measurements of each specimen' s half-wave surface 
height after testing (table I). Both specimens 
showed less favorable agreement between the data 
and the empirical solution when surface statistics 
obtained after the tests were used. 

Empirical Solution . - By comparing Figs. 5, 6, 
7, and 8 the empirical solution of Fletcher and 
Gyorog is seen to be less sensitive, like the test 
data, to differences in the specimens' half-wave 
surface heights than the theories of Refs. 1 and 2. 
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Conclusions 

Experimental data showing the variation of the 
coefficient of thermal contact conductance with 
apparent contact pressure in the light to moderate 
range (0.45x10® to 3.44x10® N/m^) are presented 
from two Armco iron cylindrical specimens having 
smooth, convex, one-half wave length contacting 
surfaces with large radii of curvature. The data 
were gathered during two loading and unloading 
sequences of each specimen at a constant interface 
temperature of 369 K and under a high vacuum condi- 
tion of 1.2>d.0" 7 torr. The data are compared with 
the theories of A. M. Clausing and B. T. Chao (j.963) 
and B. B. Mikic and W. M. Rohsenow (1966). Also, 
a comparison is made with the empirical expression 
presented by L. S. Fletcher and D. A. Gyorog in 
1970. 

The conclusions from this study are: 

1. The assumption of macroscopic elastic defor- 
mation made in both theories is inadequate for an 
accurate theoretical prediction of the variation of 
the coefficient of thermal contact conductance with 
apparent contact pressure for smooth, convex, one- 
half wave length Armco iron surfaces having large 
radii of curvature placed in contact and loaded in 
the light to moderate range. This is based on the 
cyclic data obtained and the verification that mac- 
roscopic plastic deformation of the contacting sur- 
faces occurred during the tests. 

2. Major disagreement exists between the theo- 
retical curves (based upon the surface statistics 
obtained before the tests) and the experimental 
data obtained during the first loading of the spec- 
imens. The disagreement implies that a principle 
assumption in both theories specifying the distri- 
bution of the microscopic contact spots to be uni- 
form over the contacting surfaces may be too re- 
strictive for smooth, convex, large radii of curva- 
ture specimens. 

3. Better agreement between the curves based on 
the theories and the experimental data of both 
specimens resulted when the surface statistics ob- 
tained after completion of the tests were used in 
the theoretical calculations. However, the use of 
these surface statistics is not permissible since 

a plastic set occurred during the tests. 

4. Correlation of the test data to within 35 
percent of the empirical expression of Fletcher and 
Gyorog was obtained when the surface statistics 
obtained from both specimens before the test were 
used. Like the specimens' test data, the empirical 
solution is less sensitive to differences in the 
specimens' half-wave surface heights than the solu- • 
tions based on the theories. 
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TABLE I. - ARMCO IRON SPECIMENS STATISTICS OBTAINED 
BEFORE AND AFTER CONDUCTANCE TESTS 




Upper and Lower Portion 
of Armco Iron Specimen 


2.54 diam 



Specimen 1 in Centimeters 


Before 

test 

Specimen 

d 

A. A. 

( ^max We 

Top 

46 . 2x10 "6 

2.0>CL0" 6 

35.1>0.0- 6 

Bottom 

53. 2x10 " 6 

1.4xl0‘ 6 

34.8)0.0 

After 

test 

Top 

381x10 " 6 

5. 7xl0- 6 

152. 4x10 " 6 

Bottom 

263. 3x10 “ 6 

2.5xl0- 6 

62.2x10-6 


Specimen 2 in Centimeters 


Before 

test 

Specimen 

d 

A. A. 

^ ^max^ave 

Top 

132.6X10" 6 

3.8x10-6 

64. 1x10 ' 6 

Bottom 

50.8X10" 6 

11. 6X10- 6 

99.4X10" 6 

After 

test 

Top 

282X10- 6 

2.8x10-6 

136x10-6 

Bottom 

173. 2xl0- 6 

5.1x10-6 

180 . 2x10-6 





















Figure 1. - Macroscopic constrictive resistance model of Clausing and Chao 
(ref. 1). 
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Figure 2. - Internal view of bell jar. (All dimensions are in cm. ) 
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versus apparent contact pressure for Armco Iron - 
specimen 1. Interface temperature = 369 K. 



COEFFICIENT OF THERMAL 



Figure 5. - Coefficient of thermal contact con- 
ductance versus apparent contact pressure 
for Armco Iron specimen 1. Based on surface 
statistics obtained before the test was run. 
Temperature of interface = 369 K. 
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Figure 6. -The coefficient of thermal contact 
conductance versus apparent contact pres- 
sure for Armco Iron specimen 2. Based on 
surface statistics obtained before the test 
was run. Temperature of interface = 369 K. 




Figure 7. - Coefficient of thermal contact con- 
ductance versus apparent contact pressure 
for Armco Iron specimen 1. Based on surface 
statistics obtained after the test was run. 
Temperature of interface = 369 K. 



Figure 8. - Coefficient of thermal contact 
conductance versus apparent contact 
pressure for Armco Iron specimen 2. 
Based on surface statistics obtained 
after the test was run. Temperature of 
interface = 369 K. 
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